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2 ! Abstract. We derive the requirements that a generic axion-hke field has to satisfy in order to 

^ ' play the role of the inflaton field in the warm inflation scenario. Compared to the parameter 

C^ ■ space in ordinary Natural Inflation models, we find that the parameter space in our model 

is enlarged. In particular, we avoid the problem of having an axion decay constant / that 

^^ \ relates to the Planck scale, which is instead present in the ordinary Natural Inflation models; 

in fact, our model can easily accommodate values of the axion decay constant that lie well 

CN , below the Planck scale. 
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1 Introduction 

Inflation [1-6] provides a mechanism for generating the inhomogeneities observed in the 
Cosmic Microwave Background Radiation (CMBR) [7-11], and explains the observed flatness, 
homogeneity, and the lack of relic monopoles that posed severe problems in the standard Big- 
Bang cosmology [12, 13]. Realistic microphysical models of inflation, in which the expansion 
of the universe is governed by the energy density of the inflaton field (p, are complicated by the 
requirement that the inflaton potential U{cj)) be very flat in order to explain the anisotropies 
observed in the CMBR [14]. 

Natural Inflation [15-17] is a viable model in which the inflaton is identified with an 
axion-like particle: in fact, the shift symmetry — )• + const, present in axionic theories 
assures a flat inflaton potential. Although Natural Inflation is well-motivated and it is con- 
sistent with the WMAP measurements [17], it is not an easy task to embed this model in 
fundamental theories like string theory [18], the main complication coming from the fact that 
the energy scale / at which the shift symmetry spontaneously breaks must be / > 0.6 Mpi, 
in order to agree with the constraints on the scalar spectral index n^. 

In this paper we show that the energy scale / can be as low as the Grand Unification 
Theory (GUT) scale Aqut ~ 10^^ GeV if Natural Infiation is considered in the context of 
Warm Inflation ([19]; see also Refs. [20-24]); we refer to this as the Natural Warm Inflation 
(NWI) model. 

Since axion-like particles arise in generic four-dimensional models [25, 26], string theory 
compactifications [27], and generic Kaluza-Klein theories [28], and possess attractive features 
for inflation models like a flat potential already embedded in the theory, such particles have 
been extensively discussed in the inflation literature [29-31]. In particular, an early attempt 
at lowering the Natural Inflation scale / in the Warm Inflation scenario has been discussed in 
Ref. [30] . However, the authors in Ref. [30] treat the temperature T of the primordial plasma 



as an independent variable whereas, as remarked below in this paper, the temperature T is 
not an independent quantity once the energy density of the relativistic species pr is specified 
[19]. 

This paper is organized as follows. We first fix our notation for the Warm Inflation 
scenario and for the axion particle physics in Sections 2 and 3, respectively. In Section 4 we 
analyze the dynamic of the inflaton field using the slow-roll conditions and number of e-folds 
for sufficient inflation, and we study the parameter space of our NWI model. In Section 5 
we discuss the bounds on the NWI model resulting from the measurements of cosmological 
parameters, in the light of the Wilkinson Microwave Anisotropy Probe (WMAP) mission 
plus baryon acoustic oscillations (BAO) and supernovae (SN) data [32]. Finally, discussions 
and conclusions are drawn in Section 6. 

2 The warm inflation scenario 

In the warm inflation scenario the inflaton field appreciably converts into relativistic matter 
(from here on referred to as "radiation") during the inflationary period. This mechanism is 
parametrized by the appearance of a dissipative term F in the dynamics of the inflaton field. 
In the following we reasonably assume [19] that radiation thermalizes on a time scale much 
shorter than 1/F. The energy density in radiation is thus 

p, = ^5*(T)^^ (2.1) 

where g*{T) is the number of relativistic degrees of freedom of radiation at temperature 
T. Here we do not specify g*{T) in the equations, but in the figures we will always use 
g*{T) = 228.75, corresponding to the number of relativistic degrees of freedom in the minimal 
supersymmetric Standard Model. 

Warm inflation is achieved when thermal fluctuations dominate over quantum fluctua- 
tions, or [19] 

H{T) < T, (2.2) 

where H is the Hubble expansion rate at temperature T. The effectiveness at which the 
inflaton converts into radiation is measured by the ratio 

Q = ^; (2^3) 

for Q ^ 1 a strongly dissipative regime is achieved, while Q < 1 represents the weak regime 
of warm inflation. Throughout this Section we will present the general equations for warm 
inflation, focusing on the case Q ^ 1 in the subsequent Sections when specified. 

In the following, we model the inflaton field with a scalar field (p = <i){x) minimally 
coupled to the curvature and moving in a potential U = U{(j)). The evolution of the inflaton 
field in a Friedmann-Robertson- Walker metric is described by 

+ (3/7 + F)(/. + [7^ = 0, (2.4) 

where a dot indicates the derivation with respect to the cosmic time t and U^ = dU/d(j). The 
conservation of the total energy of the system imposes that the radiation energy density pr 
satisfies 

Pr + 4Hpr = T4>'^, (2.5) 



with the term on the RHS of Eq. (2.5) describing the effectiveness of conversion of the inflaton 
field into radiation. 

The total energy density of the system at any time is ptot = 9^/'^ + U{(p,T), where 
[/((/>, T) is an effective potential that accounts for temperature effects. As discussed in 
Ref. [33], a requirement for building a consistent model of warm inflation is that finite temper- 
ature effects on the inflaton potential must be suppressed. This requirement makes possible 
to separate the effective potential of the inflaton as U{(j),T) = U{(p) + Pr{T) = U + pr- We 
thus have 

Ptoi = 1^4>'^ + U + Pr, (2.6) 

and the corresponding Friedmann equation reads 

where M-p\ = 1.221 x 10^^ GeV is the Planck mass. 

Inflation takes place when the potential U is approximately flat and the potential energy 
dominates over all other forms of energy, so that the Hubble expansion rate H is constant. 
During this period, which is known as the slow-roll regime of the inflaton field, higher deriva- 
tives in Eqs. (2.4) and (2.5) can be neglected. In formulas, 

4><^H(j), and pr <^ H pr- (2.8) 

In this regime Eq. (2.4) reads 

(/.- ^^, (2.9) 

where here and in the following we use the symbol "~" for an equality that holds only in the 
slow-roll regime. Eqs. (2.5) and (2.7) in the slow-roll regime respectively read 

Pr^^<P', (2.10) 



and 



n 8ttG 



H'^^-U, (2.11) 

so that a shallow potential U gives rise to a nearly constant expansion rate H. 

The slow-roll regime can be parametrized by a set of slow-roll parameters e, r] and /3, 
defined by 

1 ^^'A^ ,.j^%, ,.^r-^). (.1.) 



levrG \U J ' ' SvrG U ' "^ SvrG \TU ^ 

To assure the conditions expressed in Eq. (2.8), we first differentiate Eqs. (2.9), (2.10) 
and (2.11), obtaining 

He , , 



Hcj) 1 + Q V l + Q 

^'' -• ^ ^27?-/3-e + 2^^). (2.15) 



Hpr 1 + Q V 1 + Q 



Since the right hand side in the last three expressions gives the size of the terms neglected 
in the slow-roll approximation, we meet the conditions in Eq. (2.8) by demanding that these 
terms be small (see also Refs. [33-35]), 

e<l + Q, |r?|<l + Q, |/3|<1 + Q. (2.16) 

Eq. (2.16) is a generalization of the slow-roll conditions in the cool inflation that takes into 
account the parameter Q\ when Q <C 1 the dissipation term can be neglected and the slow-roll 
conditions reduce to the usual requirements in the cool inflation. 

In the warm inflation scenario, we can constraint the scale of the potential U by com- 
bining the definition of pr in Eq. (2.1) with the requirements that U ^ pr and Eq. (2.2) to 
obtain the constraint 

^'^' « f .o ^^L S ' ^Pi = 5-575*"'/'(T) X 10^8 GeV. (2.17) 

This is not a stringent bound, since in most theories of grand unification the inflaton potential 
is related to the unification scale, U^'^ ~ 10^^ GeV. 

3 Axion-like particles 

The axion [36-38] is the pseudo Nambu-Goldstone boson associated with the breaking of the 
Peccei-Quinn symmetry proposed as a solution to the strong CP problem. The theory of the 
QCD axion predicts that the axion mass at zero temperature rua is inversely proportional to 
the axion decay constant fa, with the constant of proportionality related to the QCD scale 
^QCD- In formulas, 

AqCD ~ Vl^a fa- (3.1) 

A generic axion-like theory is a generalization of the original axion theory, in which the mass 
m0 and the decay constant / of the axion theory are not related as in Eq. (3.1) but rather 
by a more generic relation, 

A = \AW, (3.2) 

where the scale A is fixed by the underlying theory. In fact, axion-like particles arise nat- 
urally whenever an approximate global symmetry is spontaneously broken, like in generic 
four-dimensional models [25, 26], string theory compactifications [27] and generic Kaluza- 
Klein theories [28]; all of these theories predict specific values for A. Axion-like particles 
mainly differ from the original QCD axion because the energy scale A differs from Aqcd^ 
depending on the underlying theory, the scale A assumes values that range up to the GUT 
scale Aqut ~ 10^^ GeV. Here the scale A is not fixed, so that the parameters m^ and / 
are unrelated [39], and Eq. (3.2) merely defines A once both m^ and / are given. 

It is well known that axion-like particles serve as suitable candidates for the inflaton 
field because the axion potential 



U = Ui<p) = A^ 



1 + cos ( A^ ^ 



(3.3) 



with A^ integer, is naturally flat whenever (p/ f is not close to ir/N. Moreover, in a generic 
inflaton model it is required that the coupling term A^ multiplying the quartic self-interaction 



0^ in the inflaton lagrangian satisfies A^ < 10 [14]. From Eq. (3.3), the quartic self- 
interaction term in the NWI model is 

2 

A* = ^. (3.4) 

In Section 5.1 we will prove that the self-interaction term in Eq. (3.4) fulfills the requirement 
from observations, see Figure 4. 

Axion-like theories provide a model for the decay rate of (j). In fact, an axion-like particle 
(j) couples to a gauge field F via the Lagrangian term 

C^F = ^(t>FF, (3.5) 

327r^ 

where F is the dual of the gauge field F and g \s a. coupling constant. For example, in the 
invisible axion theory F can be either the gluon or the electromagnetic field. In the latter 
case, axion-like particles decay in a pair of photons with a rate (see Ref. [40]) 

r = 5'y|. (3.6) 

Due to instanton effects, axion-like particles have a temperature dependent mass as [41] 

T< A, 




m^{T) = m^ { ^ ^^i ^^^ (3.7) 

In principle, this leads to a non-zero value of F^ = dT/d(j) when T > A. However, as we will 
show in Section 4.3, the NWI model is consistent only with the case T < A where the axion- 
like mass is constant. Thus, temperature effects on the axion-like mass can be neglected; in 
the following we set F^ = and, consequently, /3 = 0. 

4 Warm natural inflation 

From now on we identify the axion-like particle (p with the inflaton, using the flat potential 
in Eq. (3.3) to describe the dynamics for the inflaton field in the early universe. Here, we 
will take A^ = 1 in Eq. (3.3), so that the potential has a unique minimum at (j) = nf; the 
infiaton field relaxes towards this minimum in its evolution. In the following, we indicate the 
inflaton with (p, bearing in mind that we have assumed that (p is also a pseudo-scalar particle 
for which the formulas in Section 3 apply. 

Furthermore, from here on we are interested only in the strongly dissipative regime of 
warm inflation, Q ^ 1. In this limit, the only parameters in the theory are the inflaton mass 
m^, the decay constant / and the dissipation term F. In table 1 we have collected the most 
important quantities of the theory in terms of these parameters 

We now examine the constraints on the NWI model coming from the slow-roll conditions 
in Eq. (2.16), the requirement for sufficient inflation, and the WMAP measurements on the 
power spectrum of density and tensor perturbations. 



Quantity 



Equation(s) used 



C/((/))=m2/2(l + cos<^//), 
U^((t>) = - mlf siiKp/f, 

™^/ /^(l + cos0//), 



H 






m^/_\/ 247r(l+cos 0//) ' 
P*^ - V 128^ r ^i+cos<^//' 



32^ rg^T) yi+cos «!>//' 

_ 1 sin^0// 

^ ~ 167r G /2 (l+cos(/>//)2' 

^^ 1 cos </)// 

^ ^ ~87rG/2 1+cos (/>//' 



Eq. (3.3); 

derived from Eq. (3.3); 

Eqs. (2.11) and (3.3); 

Eqs. (2.9), (2.11) and (3.3); 
Eqs. (3.3) and (4.1); 

Eqs. (2.9), (2.10), (2.11) and (3.3); 

Eqs. (2.1), (2.9), (2.10), (2.11) and (3.3); 

Eqs. (2.12) and (3.3); 
Eqs. (2.12) and (3.3). 



Table 1. Expressions for sonic derived quantities in the theory, vahd during slow-roU and Q ^ 1. 

4.1 Slovif-roll conditions 

During warm inflation the values of the slow-roll parameters e and r] must both be smaller 
than Q, see Eq. (2.16) with Q ^ 1. Here we neglect the condition /3 <C Q because in this 
model /3 = 0. Inflation ends when one of these two slow-roll conditions is violated. 
Writing the parameter Q during slow-roll as 



Q — - 

V247rG[/' 

the slow-roll conditions e <^ Q and t] <^ Q read 



1 m, 



< 



WttG \U J V247rGC/ 



and 



7] 



Urh^ r 



(4.1) 



(4.2) 



(4.3) 



SttGU V247rGC/ 

Using the expression for the potential U{(j)) in Eq. (3.3), these slow-roll conditions give 
respectively 

1 - cos 0// 



and 



y/l+COS(j)/f 
- COS 0// 



^ a, 



< 



a 



^1 + cos cP/f 2 
where we have defined the combination 



_ /327r Tf ,^,j. f -1 



V9' 



(4.4) 



(4.5) 



(4.6) 



Here and in the following, a number y indexing some quantity with units of energy indicates 
that such quantity has been divided by 10^ GeV: for example, /ig = //lO^^GeV. 

The slow-roll regime ends when the field (p reaches a value (j)f for which one of the 
conditions e <^ Q or ij <^ Q is no longer satisfied. We have checked numerically that in 
general the second condition r] <^ Q is tighter than the first one e <^ Q, but for a ^ 1 the 
two conditions give the same result, 



/ u 



V8 



a 



(4.7) 



Since observations favor a large value of a > 10 (see below), in this paper we use Eq. (4.7) 
to define (pf. 

4.2 Number of E-folds 

The number of e-folds is defined as 



/■t2 

Ne = ln(a2/ai) = / Hdt, 
Jti 



(4.8) 



where ai and 02 are the values of the scale factor a(t) appearing in the Friedmann metric 
when inflation begins and ends, respectively. Sufficient inflation requires 



N, > 60. 



(4.9) 



During the inflationary stage, the value of the inflaton field decreases from the initial value 
(pi to the value at the end of inflation (pj defined via Eq. (4.7). We now derive a relation 
between (pi and (pf, using the definition of A'^e above. In the case of a slow-rolling of the 
inflaton and in the strongly dissipative regime, Eq. (4.8) reads 



iV. 






(4.10) 



where we used Eq. (2.9) with T ^ 'iH. Using the axion-like potential in Eq. (3.3) and the 
Friedmann equation we obtain 



[8^ Tf f't'f/f , Vl + cosx a , X 

Ne = \ T-;— / dx = —p= m tan — 

V 3 m^Mpi J^jf smx V2 4 



<^>^/f 



or 



tan 



4/ 



tan -— Exp 
4/ 



V2Ne 



a 



In general, we use Eqs. (4.7) and (4.12) to obtain the value of (pi/ f, 



(4.11) 



(4.12) 



/ 



4 arctan 



1 — tan 



1 -|- tan 



f^Exp 

V2a 



V2N, 



a 



(4.13) 



The appearance of trigonometric functions is due to the shape of the potential U{(p), that 
differs from a pure quadratic one and contains a cosine function itself. 
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Figure 1. The initial value of the angle 0i// as a function of a = •\/327r/3 Tf/m^Mpi for different 
value of the numbers of e- folds N^. Blue dotted: Ne ~ 40; Red dashed: N^ = 60; Green dot-dashed: 
A^e = 80. Also shown is the line (j>i = tt/ (Blaek dashed line). 



In Figure 1 we show the value of (pi/ f given in Eq. (4.13) as a function of a. Considering 
the three parameters m^, / and T from which a depends, moving towards greater values of 
a corresponds to a decrease in m^ or to an increase in / or F, once the other two parameters 
have been fixed. 

For larger values of a, the initial angle (pi/ f approaches vr and the axion-like potential 
is not distinguishable from a pure quadratic one. This fact is shown in Figure 2, where we 
compare the value of the angle (pi/f with A'^e = 60, together with the value of (pi if the inflaton 
potential were a pure quadratic one. 



A^ 



instead of Eq. (3.3). For a > N^, it is 



vr - NeVs/a, 



(4.14) 



(4.15) 



and the dynamics of the field can no longer discern between the two potentials U{(p) and 

Uqnad{(p)- 

4.3 Parameter space of the NWI 

We have derived our results in the strongly dissipating regime F ^ 3H: this constraints the 
parameters as 



a S> 167r 



/ 



Mpi 



) yrr 



COS( 



(4.16) 



This constraint is not severe for values of / around the GUT scale, which is the region of 
interest in this paper. In fact, for / ~ 10^^ GeV a numerical solution of Eq. (4.16) is 



a > 4.8 X 10~^ /fg. 



(4.17) 
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Figure 2. The initial value of the angle (t>i/ f as a function of a = ■\/327r/3 T J / m^Mp\ for the case of 
the axion-like potential in Eq. (3.3) (black dotted line) and for the quadratic potential in Eq. (4.14) 
that approximates Eq. (3.3) around 0^ ~ tt/ (black dot-dashed line). We fixed N^ = 60. 



independently of Nf^. This bound is much more loose than the requirement of flatness of U{(j)) , 
which yielded a > 1. When / ~ Mpi, the constraint in Eq. (4.17) starts depending mildly 
on the number of e-folds: for example, taking / = Mpi in Eq. (4.16) gives a ^ 20A''g'^^. 

We now discuss the condition T < A, that allows us not to consider temperature effects 
in the axion mass in Eq. (3.7). Using the equations in Table I for T and A we obtain 



15 sin^(/)i// 
a > — :t 



vr^ 



a/1 + cos (pi/f 



(4.18) 



which is always satisfied for all values of a, as can be checked numerically. Modeling T 
through Eq. (3.6) gives E^ = and the slow-roll parameter /3 = in this model. Summing 
up, the temperature of the plasma in the NWI model is constrained hy H < T < A, with 
the additional constraint F ^ 3H that yields to Eq. (4.16). 



5 Perturbations from inflation 

As mentioned in the Introduction, one attractive feature of inflation is that scalar and tensor 
perturbations emerge during this epoch. These features later evolve into primordial fluctua- 
tions of the density profile and gravitational waves, that might leave an imprint in the CMBR 
anisotropy and on the large scale structures [7-11]. Each fluctuation is characterized by a 
power spectrum and a spectral index, respectively A^(fc), n^ for density perturbations and 
A^(fc), riT for tensor perturbations. Here we use results from the warm inflation literature 
to set constraints on the NWI model. 

5.1 Scalar power spectrum 

The spectrum of the adiabatic density perturbations generated by inflation is specified by the 
power spectrum A^(/c), which depends mildly on the co-moving wavenumber k according to 



a spectral index ns{k) as [42] 



2 „., _ k' Pitik) _ ,2„> fky""-' 



^n(k) ^ -^^ = A4(».„) ^jjj . (5.1) 

The function A^(A;) describes the contribution to the total variance of primordial curvature 
perturbations Tt due to 7^ at a given scale per logarithmic interval in k [32]: the WMAP 
collaboration reports the combined measurement from WMAP+BAO+SN of A|j(/i:o) at the 
reference wavenumber k = ko = 0.002 Mpc~^, 

A|(A:o) = (2.445 ± 0.096) x 10^^ (5.2) 

where the uncertainty refers to a 68% likelihood interval. The RHS of Eq. (5.1) is evaluated 
when a given co-moving wavelength crosses outside the Hubble radius during inflation, and 
the LHS when the same wavelength re-enters the horizon. In Eq. (5.1) we have used the 
notation in Ref. [32] for the density perturbations. Other authors use the symbol P-ji^k) for 
our A^(fc) and 7^| for our PTi{k), and might differ for factors of 27r^. 

In both warm and cool inflation models, the scalar power spectrum has the form 

^Uk)=i^{H)] , (5.3) 



where (j) ~ — [7^/(3/7 + T) and {6(j)) describes the spectrum of fluctuations in the inflaton 
field. The LHS of Eq. (5.3) is computed at the time at which the largest density perturbations 
on observable scales are produced, corresponding Nf, e-foldings before the end of inflation. 
Quantum fluctuations predict [8] 

rr 

('^</')quantum = TT i (5-4) 

zvr 
while thermal fluctuations provide (we used Eqs. (21) and (23) in Ref. [43]) 

(<5</')thcrmal = f .^^^3 j • (5.5) 

Using Eqs. (5.3) and (5.5) for the variance of fluctuations we obtain the scalar power spectrum 
in the strongly dissipative regime of warm inflation as 

which is the same result as in Refs. [33, 34] with an extra factor l/(27r)^ included to account 
for our normalization of the power spectrum. For comparison, the power spectrum in the 
usual cool inflation is found as 

ALooi(M=f#^l , (5.7) 



1/4 



V27r(/) 
so that density perturbations are larger in warm inflation by a factor 



Aieooi(^^o) \^ H-' J \ml,A, 



10^ -^^ • (5-^ 



10 



From now on, scalar perturbations are considered only in the warm inflation scenario, 
so we suppress the index "warm" in Eq. (5.6). Using the expressions for H and T in Table I 
we write Eq. (5.6) as 



.Ml 



sin (f)i/f 



In deriving this last expression we have used the fact that the largest density perturbations are 
produced when (j) = (j)i [15]. Another form of Eq. (5.9) that is more useful for computations 
is 

9*{T) 



A^{ko) ~ 3.7 X 10~^^ 



228.75 



V'^' 3/4^/2 Al+C0S.^.//)3/A 

J " ^12 1^ sin<^.// ) 



3/2 



(5.10) 



Eq. (5.10) defines the power spectrum in terms of a, T and N^, the latter appearing implicitly 
in the definition of 4>i. 

We equate the expression for the power spectrum in Eq. (5.10) to the measured value 
from WMAP in Eq. (5.2) to obtain a relation between T and a, see Figure 3. In Figure 5.2, 
we only consider the case A'^e = 60, while other values of N^ do not modify the curves 
sensibly unless a 3> 1. The thickness of the line approximately corresponds to the one-sigma 
uncertainty in Eq. (5.2). From Figure 3, the dissipation F reaches a maximum value Fmax at 
large values of a, which we find analytically from equating Eqs. (5.2) and (5.10) in the limit 
a > 1, 
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Figure 3. The dissipation term F as a function of a, given by Eq. (5.2) witli A^(fco) = (2.445 ± 
0.096) X 10^^. We used Eq. (5.10) for the analytic expression of the scalar power spectrum. 



To study the strength of the quartic self-interaction in the NWI model, we use Eqs. (3.4) 



and (4.6) to eliminate F in Eq. (5.10) and obtain a relation between A,^ and a, 

3/2 



A^{ko) ~ 6.89 X 10 



^14 I g*(^) V'^' 9/4x3/4 Al+COS./>,//)3/4 ' 

.228.75; ^ [ sin (/>,// 



(5.12) 



11 



The dependence of the quartic self-interaction A,^ on a thus obtained is shown in Figure 4. 
We see from this figure that the self-interaction is A^ < 10^^*^ for any a, and the NWI model 
can easily satisfy the constraint for the quartic self-interaction term A^ < 10~^ obtained in 
Ref. [14]. 
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Figure 4. The inflaton quartic self-coupling A^ as a function of the parameter a, Eq. (5.12). The 
parameters A^(fco) and N^ are fixed as indicated in the figure. 



5.2 Scalar spectral index 

The scalar spectral index n^ describes the mild dependence of the scalar power spectrum on 

the wavenumber A;, as in Eq. (5.1). We expand the spectral index around the reference scale 

fco as 

1 h 

(5.13) 



nAk') = n, H — r In — . 



where n^ = ns(/co) and the spectral tilt r is 

dugik) 



d\nk/k(j 
Using Eq. (5.1), the scalar spectral index is [42] 



k=ko 



d 



n. 



In 



A|(fc) 



dlnk/ko "^ A^{ko)' 
In warm inflation, the spectral index is ([35]; see also Ref. [44]) 



n.o 



1 



1 
Q 



;v 



r/5 



(5.14) 



(5.15) 



(5.16) 



while the spectral tilt depends on higher orders in the slow-roll parameters, and for this 
reason it will be neglected here. With the expressions in Eq. (2.12) for e and t], and with 
/3 = 0, we find 



n.. 



1 



3 2 

167rGC/2Q V^'^'^^"4^' 



(5.17) 
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or, using t/((/)) in Eq. (3.3) and its derivatives at 
rj in Table 1, 



bi, together with the values of Q, e and 



ric 



_ 3 /_3_ m^ Mpi 3 + cos (pi/f 



8V27r Tf ^1 + cos 4>,/f 



1 



1.50 3 + cos (pi/ f 



a y/l + cos 4>i/f 



(5.18) 



This expression differs from what found in the "cool" Natural Inflation scenario (see Eq. (11) 
in Ref. [17]). 

Figure 5 shows the dependence of the spectral index on a as in Eq. (5.18), for different 
values of the number of e- folds N^. We have compared the values of n^ with the combined 
WMAP +SN+BAO measurement [32], 



Tic 



0.960 ± 0.013 at 68% C.L., 



(5.19) 



so that the shaded yellow region in Figure 5 corresponds to the 68% C.L., while the light blue 
region is the corresponding 95% C.L. region. We see that small values of Nf. ~ 50 and large 
values of a > O(IOO) are favored, although values of A'^e = 60 — 70 can be accommodated in 
the 68% C.L. region for 100 ^ a < 300 or in the 95%C.L. region for all values of a. 
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Figure 5. The scalar spectral index Ug, Eq. (5.18), for different values of the number of e-folds. Blue 
dotted: N^ = 40; Red dashed: N^ = 60; Brown dot-dashed: N^ = 80. Also shown are the C.L. regions 
from the combined WMAP+BAO+SN measurement, Eq. (5.19). Yellow: 68%; violet: 95%. 



For large values of a the scalar spectral index Ug is essentially independent of a, while 
for a < A^e it is the dependence on A^e that vanishes. This behavior can be summed up by 
considering Eq. (5.18) in these two limits. 



1 



n. 



4.2/a a<Ne 
1.5/Ne a>Ne 



(5.20) 



This type of dependence of Ug on either A^e or a has been also noticed in models of Natural 
Inflation set in the usual cool inflation scenario, see for example Eq. (12) in Ref. [17]. In 
the case in Ref. [17], in which there is no dissipation term F, and the inflaton mass m^ is 
adjusted so that A ~ 10^^ GeV, the result analogous to Eq. (5.20) is expressed in terms of / 
instead of a. 
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5.3 Tensor power spectrum 

Since warm inflation considers thermal fluctuations instead of quantum fluctuations to gen- 
erate scalar perturbations, it is only density fluctuations that modify in this scenario while 
tensor perturbations show the same spectrum as in the usual cool inflation [33]. Defining the 
tensor power spectrum as 

AfW-^ = A^(M(^)"\ (5.21) 

where the tensor spectral index ut is assumed to be independent of k, because current 
measurement cannot constraint its scale dependence. WMAP does not constraint A^(fco) 
directly, but rather the tensor-to-scalar ratio 

which qualitatively measures the amplitude of gravitational waves per density fluctuations. 
The WMAP-|-BAO-|-SN measurement constraints the tensor-to-scalar ratio as [32] 

r < 0.22 at 95% C.L.. (5.23) 

Since in warm inflation the scalar power spectrum is enhanced by the quantity ^ in 
Eq. (5.8) with respect to the value in cool inflation, the tensor-to-scalar ratio and thus 
gravitational waves are reduced in the NWI model by the same amount ^. Using Eq. (5.6) 
for the scalar power spectrum and the expression for the tensor power spectrum, 

^riko) = ^, (5.24) 



the tensor-to-scalar ratio in warm inflation is 

128^ 4 



Ml, JhTT' 



'PI 
With the values in Table 1 we find 



(5.25) 



128 ^^V a'/'m L^] ^uuklf)^_ 
^'HlSoj ^* ^''[r^MlJ (l + cos,^.//)V8- 



In the limit a ^ A'^e, with (/)j ~ vr — \'sN(,/a, we obtain 

..«3.10-A.f*(|EL)"*(£ii-)\ fo.„»A... (5.27) 

Using the WMAP measure of the tensor-to-scalar ratio in Eq. (5.25), we constraint the 
inflaton mass in the NWI model to 

-<--»'^-(^)-(t)--(|ii)-, ..2S, 



14 



where we used the expression for Fmax in Eq. (5.11). Future measurements wih not substan- 
tiahy improve the bound in Eq. (5.28), because of the power 1/4 raising r. As an example, 
the forecast PLANCK measurement will constraint the tensor-to-scalar ratio by one order of 
magnitude with respect to WMAP, r < 0.01 [45], thus the bound in Eq. (5.28) will approxi- 
mately lower by a factor of two. 

Results are summarized in Figures 6 and 7. In Figure 6 we show the constraints on the 
NWI model from the WMAP data in the r-Ug plane, in the case / = 10^^ GeV. The WMAP 
data are better in agreement with a low number of e-folds N^ = 40 — 50, although values up 
to Ng = 80 can be accommodated within the 95% C.L. region. We conclude that a value of 
the axion decay constant of the order of the GUT scale / ~ Aqut can be easily embedded 
in the NWI model, making it possible to construct microscopic theories of warm inflation 
with axion-like particles at the GUT scale. For this value of the axion decay constant, the 
expected value of r and thus the amount of gravitational waves that are produced with this 
type of inflation is extremely low. If gravitational waves are found with r ~ 10"^^ or above, 
this model has to be abandoned and a Aqut valued axion decay constant / is no longer 
viable. 
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Figure 6. Prediction from the NWI model and WMAP constraints in the r-Us plane for / = 10^^ 
GeV. Blue dotted line: iVe = 40; Red dashed line: N^ = 60; Brown dot-dashed line: N^ = 80. The 
yellow and violet regions are the parameter spaces allowed by WMAP+BAO+SN at 68% and 95% 
C.L. respectively. 



For comparison, we show the constraints on the NWI model from the WMAP data in the 
r-Ug plane when / = 10^^ GeV in Figure 7. We see that the expected amount of gravitational 
waves has raised with respect to the case in Figure 6 by ten orders of magnitude, while data 
still favor relatively low values of the number of e-folds. For a given value of A'^e, the curves 
in Figures 6 and 7 rigidly move along the vertical r axis when changing the parameter /„. 

6 Discussion and conclusion 

The parameter space of the inflaton field in the NWI model is broader than that in the usual 
Natural Inflation, because of the presence of the extra quantity F describing dissipation and 
of the enlarged slow-roll conditions in Eq. (2.16). 
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Figure 7. Prediction from tlie NWl model and WMAP constraints in the r-Us plane for 
/ = IQis GeV. Blue dotted line: TVe = 40; Red dashed line: TVe = 60; Brown dot-dashed hne: 
A^e = 80. The yellow and violet regions are the parameter spaces allowed by WMAP+BAO+SN at 
68% and 95% C.L. respectively. 



The NWI model allows to lower the value of the axion decay constant / from the Planck 
scale resulting in Natural Inflation to the GUT scale; a ratio f/Mpi ~ 10^"^ helps overcoming 
some difficulties encountered in Natural Inflation model-building. 

A decay constant of the order of the Planck scale is still possible, since / is not bound 
from above. In the case / ~ Mpi, Natural Inflation and NWI can be distinguished observa- 
tionally by a measurement of the tensor-to-scalar ratio r, see Figures 6 and 7. This difference 
in the value of the ratio r comes from the fact that in the NWI model the amplitude of 
gravitational waves is suppressed by the factor ^ in Eq (5.8) with respect to the standard 
cool inflation. 

Measurements of the scalar spectral index n^ favor larger values of a > A^ei as shown 
in Figure 5. In this region, the value of (j)i in Eq. (4.13) can be approximated with (pi ~ 
vr — \/SNf,/a, see Eq. (4.15), and the axion-like potential in Eq. (3.3) cannot be distinguished 
from a pure quadratic one through the dynamic of the inflaton. In this condition, one would 
need to obtain the value of the self-interaction A^ independently, for example from considering 
density fluctuations in the CMBR; Natural Inflation and NWI models predict a precise ratio 
between the height of the potential and the self-interaction coupling constant, U/X^ = 12 /^, 
whereas other inflaton models predict different values of this ratio. 

The possibility that the axion-like energy scale / is of the order of the GUT scale when 
warm inflation is considered has been already taken into account in Ref. [30]. However, we 
differ from Ref. [30] in various points: for example, we do not consider the temperature 
as an independent variable, because in warm inflation the radiation bath is thermalized and 
temperature is expressed is linked to the radiation energy density pr as in Eq. (2.1). Moreover, 
here we have included a detailed analysis of the CMBR observables r, A^(A;) and n^ that 
was missing in Ref. [30]. 

We have presented a model in which Natural Inflation takes place within the warm 
inflation scenario: in such model, the decay constant / no longer ties to the Planck scale as 
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in the usual Natural Inflation model, but it can be as low as the GUT scale, / ~ 10^^ GeV. 
We have shown the viability of the NWI model and its agreement with current astrophysical 
data. 

During the completion of the present paper, the work in Ref. [46] where a similar model 
is considered came to our attention. The authors in Ref. [46] constrained the Natural Warm 
Inflation model using data from the WMAP experiment and obtained values of / that are 
generally below the Planck scale, in agreement with the results in this paper. 
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